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en ■ Abstract 

^ ' Fractal structures and non-Gaussian velocity distributions are characteristic prop- 

Cn , erties commonly observed in virialized self-gravitating systems such as galaxies and 

interstellar molecular clouds. We study the origin of these properties using a one- 
dimensional ring model which we newly propose in this paper. In this simple model, A'' 
particles are moving, on a circular ring fixed in three dimensional space, with mutual 
interaction of gravity. This model is suitable for the accurate symplectic integration 
^^ ] method by which we argue the phase transition in this system. Especially, in be- 

tween the extended-phase and the coUapsed-phase, we find an interesting phase {halo 
-phase) which has negative specific heat at the intermediate energy scale. Moreover 
in this phase, there appear scaling properties and non-thermal and non-Gaussian ve- 
Q ' locity distributions. In contrast, these peculiar properties are never observed in other 

5— ( ' gas and core phases. Particles in each phase have a typical time scale of motions 

C/3 , determined by the cutoff length ^, the ring radius R, and the total energy E. Thus 

Cu ' all relaxation patterns of the system are determined by these three time scales. 
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I. INTRODUCTION 

Many astrophysical objects in our universe consist of mutually interacting elements through gravity. 
If they are almost isolated systems, they are called self-gravitating system (SGS). For example, galaxies, 
clusters of galaxies, globular clusters, and molecular clouds are thought to be typical SGSs. Their 
statistical properties are often characterized by non-Gaussian velocity (or pair-wise velocity) distributions 
[Ql , fractal structures |0] , and the scaling relation between the mass density and the system size [g| . Most 
of these objects are thought to be gravitationally virialized. Therefore pure gravitational force seems to 
play an essential role in characterizing the above statistical properties of SGS independently of initial 
conditions. 

There have been some theoretical approaches to explain the fractal structures in SGS from the viewpoint 
of criticality and phase transition in gravo-thermo dynamics M. Strictly speaking, the ordinary SGSs 
in three-dimension (3-D) cannot attain genuine stable equilibrium, because the gravitational force does 
not vanish at long distances (IR-divergence) and diverges at short distances (UV-divergence). These 
properties of gravity cause gravothermal catastrophe in a self-gravitating gas system enclosed by a solid 
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adiabatic wall. In fact, the isothermal sphere is not always stable, since the entropy does not necessarily 
take the local maximum for this configuration |^-0 . Therefore the introduction of a small-scale cutoff as 
well as a large-scale cutoff is inevitable for such unstable systems to attain the final equilibrium. Though 
the introduction of the cutoff prevents the gravothermal catastrophe, 3-D gravitational system has a 
phase with negative specific heat when treated in the microcanonical ensemble pi . Even in the canonical 
ensemble, this system is highly unstable. When the temperature decreases, the system experiences 
the violent first-order phase transition from the gas phase into the cluster phase. According to these 
arguments, no stable equilibrium states are theoretically expected in the 3-D gravitational system with 
or without cut-off in any ensemble. 

On the other hand, in the real world, we do not need the stable equilibrium states for describing 
SGS; meta-stable equilibrium states do actually appear in the dynamical description such as coUisionless 
Boltzmann equation. The finite lifetime of such structures is sufficient to explain the present structures 
of SGS, even if they are expected to evolve further into different quasi-stable states through two-body 
relaxation. 

The relaxation process of SGS has been mainly discussed in the one-dimensional gravitational sheet 
model (IDS) pO| , pl| . In this model, many parallel sheets interact with each other through constant force 
which never decays at distant places. Though the interaction is long-ranged, no phase transition occurs 
in IDS. Thermodynamics of IDS is exactly solved and actually, in numerical calculations, the system 
reaches thermal equilibrium long after it attains virial equilibrium. In this model, the virial condition 
gives the relation 2 {K) — {V) between the time averaged kinetic energy {K) and the potential energy 
{V). Therefore, contrary to the 3-D SGS, specific heat of IDS is always positive. Thus the relaxation 
process in IDS would be quite different from that in 3-D SGS. Another well-known one-dimensional 
model which has long-range force is the Hamilton Mean Field (HMF) model, in which phase transition 
does occur ||l^. There have been a lot of studies on the relaxation process of HMF. Actually, Levy type 
jumping motion of constituent particles and super-diffusion process have been revealed in HMF |l3]. In 
this HMF, though the specific heat becomes negative in the quasi-stable state 0, it remains positive in 
the thermal equilibrium state and the phase transition turns out to be the second-order. Despite some 
common properties of the models, the interaction form in HMF is also quite different from that of real 
3-D gravitational systems. 

Several other simple models have been proposed |1^] in order to characterize the 3-D SGS much faith- 
fully. For example, the cell model and its extended version are the simplest models which show phase 
transition. In these models, the pair interaction potential changes its value only at the cell boundary 
||16| , |l7| . At low temperatures, most particles are trapped within several clusters, and at higher temper- 
atures, these clusters melt and the particles can move freely. However, it seems difficult to examine 
the relaxation process of this system numerically, since the interaction is not efficient to cause sufficient 
relaxation; particles have no interaction within each cell. There are also several numerical analysis with 
the different type interactions including the lower and higher cutoffs |l^Jl^. However, the value of the 
lower cutoff used in realistic numerical calculations of these models seems to be too large to extract the 
intrinsic properties of realistic 3-D SGS. 

Thus it seems indispensable for us to have a model which faithfully refiects the characteristic properties 
of 3-D SGSs; the model has negative specific heat and shows the phase transition similar to 3-D SGSs, and 
is capable for us to analyze quasi-equilibrium states which would be realized before the system reaches 
the complete equilibrium. 

In this paper, therefore, we first propose a new model (Self Gravitating Ring: SGR) which a) has neg- 
ative specific heat, b) shows phase transition representing 3-D SGS, and c) is numerically tractable. This 
model consists of iV-particles, which can freely move on a ring with a fixed radius, mutually interacting 
through 3-D gravity. An excellent point of this model is that the force is genuine 3-D gravity while the 
calculation is essentially one-dimensional. Moreover, the Hamiltonian permits the accurate symplectic 
integration method by which we can analyze the nature of quasi-equilibrium states and phase transitions 
within very long time scales in this model. 

In section O, we introduce the SGR model, and in section [II identify three quasi-equilibrium states 



including the state with negative specific heat. In section IV, we analyze the particle motions from a sta- 
tistical point of view, and then, we study the relaxation process toward the thermodynamical equilibrium 
in section 0, and the scaling structure in section VI . Finally, we discuss our results in section |VII| . 



II. SELF-GRAVITATING RING MODEL 

In this section, we introduce the self-gravitating ring (SGR) model, in which particles interact with each 
other through genuine 3-D gravity while the particle motion is constrained on a 1-D ring. By utilizing 
this model, we study the gravitational phase transition, which was not feasible in the IDS model |ll[] . 
The lack of phase transition in IDS is related with the fact that the potential of IDS increases linearly 
without bound, and therefore there is no characteristic energy scale required for a particle to escape from 
the cluster. On the other hand, we have the characteristic energy scale to bind particles in our SGR 
model. 

We consider a system of self-gravitating N particles with mass to, whose motions are smoothly con- 
strained on a circular ring with a fixed radius R without friction (Fig.nl) . Each pair of particles interact 
with each other through 3-D gravitational force. The distance between a pair is measured by the length 
of the straight line combining the pair and not by the minor arc of the ring. 




FIG. 1. The SGR model with a fixed radius R. The particle locations are specified by the angles measured 
from a fixed direction. A pair of particles at 9i and at 6j interact with each other through the inverse-square-law 
3-D gravitational force; the distance is measured by the straight line in the picture and not by R\6i — Oj\. 

The Hamiltonian of this system becomes 

P 2TOi?2 ^ ^ ^ V2i?Vl-cos(0.-0,) + 6- ^ ■ ' 

The position of z-th particle is fully described by the angular variable 6i asr,i — {R sin 0i,R cos 6i). The 
momentum conjugate to 0i is given by Pi — niR^dOi/dt. The UV-cutoff parameter e truncates the 
diverging gravitational force at around the distance ^ = ^/2eR. 

We first introduce three dynamical time-scales which are apparent in the above Hamiltonian. They 
are parameterized by the ring radius i?, the cutoff scale ^, and the total energy E of the system. 

When the system is almost uniformly filled by moving particles, we have the longest dynamical time 
tfi defined by 



i^=W-^. (2.2) 

^ V GNm ^ ' 

During this time t/j, a typical particle goes around the ring once. Therefore, the time-scale for the whole 
system to attain thermodynamical equilibrium, if any, is at least larger than this time scale. 

When all the particles collapse completely into a core, we find the shortest dynamical time-scale i^ 
defined as 



*^^^/gL- (^-^^ 



During this time t^, a typical particle bound to the core oscillates once, n 

There is an another time-scale in between the above two extreme time scales. Suppose the system is 
stably confined in a region with a scale r which satisfies 

^«r«R. (2.4) 

For this scale r, the ring can be approximated as an infinite straight line and the cutoff ^ can be neglected. 



Then the leading term of the denominator of the potential term in the Hamiltonian (2.1) becomes \9i — 9j\ 
and 
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In this case, the potential term Vp satisfies Euler's theorem for homogeneous function, 

E^#=-^P' (2-6) 

and the ordinary virial condition holds 



.1 ^^» 



2 (Kp) = - (Vp) , (2.7) 

where angle brackets represent the long-time average. [| According to this virial relation, the typical size 
of a system r is related with the total energy E as 



Hence the above condition for r (2.4) reads 



^T >> "" >> ^^- ^'-'^ 

Moreover from the virial relation, the velocity dispersion is given as 

2(i^p) = -(yp)=m7V(w2), (2.10) 

where v = dr/dt and therefore 



VW)-J^. (2.11) 



iN 



From this equation (2.11), the crossing time is defined as 



r Gm^l^N^I^ 

tE = , = 7= ^- (2.12) 

ypy 4V2|£;|3/2 ^ ' 

This is the intrinsic time scale associated with genuine gravity independent of the cutoff ^ and of the 
system size R. [] 



* In the next section, we classify all the particles into three species, core , halo , and gas . The typical time 
scales tR and t^ respectively characterize the gas and the core species. 

^ Outside of the above region in SGR model, this form of virial relation would be modified even in the limit of 
£, -^ 0,as is shown in Appendix A. 

* As we will show in SecffM, this is also the time scale for the particles in halo species in SGR model. 



The above introduced cutoff parameter ^ connects the two limiting cases in the following sense. In the 
limit ^ ^ 0, the SGR model becomes genuine 3-D gravity at small scales, while in the limit £, ^ R, the 
model almost becomes the HMF model. The latter is because, in this limit, the shape of the potential 
around 9 — kn{k — 0, ±1, • • •) becomes almost identical as does in the HMF model ^. 

For numerical simulations, we need to make all physical variables non-dimensional; we use m, R, and 



tii for the unit of mass, distance, and time, respectively. In these units, the physical Hamiltonian (2.1) 
reads 



Hp = Gm^—H, (2.13) 

R 



where 



N N 

H^IYpI-Y— ^ ^ (2.14) 



^/2N^ - cos 



2 — 1 1<J ^ V *J 

The dimensionless momentum pi is given by pi = dOi/dr and the dimensionless time r is introduced as 
T = t/tn. The above form of Hamiltonian permits us to use a powerful symplectic integrator p9| , with 
which the total energy is conserved with extremely high accuracy, even beyond thousands of dynamical 
time. Typical magnitudes of errors for the total Hamiltonian H{t) and the total momentum P{t) in 
our simulations up to r - lO'^, are {H{t) - H{0))/H{0) - 0(10"^) and P{T)/Pr.m.s{T) ^ 0(10"^), 

respectively, where Pr.m.sir) = (j^iLiPi) 

III. CLASSIFICATION OF PHASES AND PARTICLES 

We now study the quasi-equilibrium state that appears in a transient stationary stage in our SGR 
system. Though this state is not absolutely stable, it generally appears in SGSs during sufficiently long 
time before the system finally approaches the equilibrium state characterized by the equipartition of 
particle energy. 

We would like to extract universal properties observed in this transient state; only the transient de- 
scription is possible and necessary to explain observations. 

A. Negative specific heat of SGS 

First, we study the phase diagram of SGR model in the temperature-energy plane. The T-U relation 
is shown in Fig.0 and we observe that the region with negative specific heat, i.e. negative slope region, 
apparently exists. The temperature T and the internal energy per particle U of the present system are 
defined respectively by 

t/.f . (12) 

The internal energy U is conserved in the present microcanonical system, while the temperature T is 
defined as the twice of the time-averaged (represented by angle brackets) kinetic energy per particle. The 
phase diagram in FigJ3 is time averaged until t — 2 where the system archives virial equilibrium (see 
Fig.y). Throughout this paper, except explicitly mentioned, we fix the total number of particles N = 100 
for simplicity and therefore the remaining relevant parameters which characterize the quasi-equilibrium 
state of the system would be e and U. 



The exact HMF model is reproduced in the limit ^ ^ oo. 
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FIG. 2. Diagram of temperature T vs. energy per particle U for three different cutoffs e — 2.5 x 10 ^, 
1.0 X 10~®, and 2.5 x lO"'^. In each T-U curve, there are two critical energy scales; t^c-low(^) ^^'^ '-'^c-hiffh^^)' 
between which negative specific heat dU/dT < appears. 



From Fig.g, we find two characteristic energy scales where dT/dU — 0; U^,_■[Q^{e) at a low energy side 
and t^c_high(^) ^^ ^ high energy side. The energy scale U^_-^^„^{e) corresponds to a mean gravitational 
binding energy per particle, which is estimated as {l/{y/2N\/l -cos6' + e)) x {N{N - l)/2)/N ~ 0(1), 
while t/r._]r,^,r(e) strongly depends on the cutoff e. This cutoff dependence can be estimated from the 
condition (^.9|) which, in our normalization, becomes 



4V^ 



<[/< 



(3.3) 



Under this condition, the negative specific heat condition in virialized state would b e ju stified. Actually, 
substitution of three different cutoffs e = 2.5 x 10""^, 1.0 x 10~^, and 2.5 x 10~^ into (3^) yields the lower 



limits of (3.3) 353, 177, and 35.3, respectively. Thus the condition (3.3) correctly describes the region of 
negative specific heat in Fig.H with sufficient accuracy. However, the slope of each T-U line in Fig.g is 
less steep than the value —2 which is expected for the virial c ond ition of 3-D gravity. This discrepancy is 
probably because the energy range satisfying the condition (3.3) is too narrow for the ideal T-U curves 
according to our choice of the cutoff parameters; even the smallest cutoff we took may not be enough to 
make the relation converge. 

In the range between these two energy scales, i.e. C^c-low(^) < U < C^c_high(^)' ^^'^ system has 
negative specific heat, which suggests the existence of phase transition ||q,[l6|-[l^ . Actually in the system 
with negative specific heat, a slight energy injection from outside decreases the system's temperature and 
induces further energy flow from outside. Then this catastrophic temperature reduction induces rapid 
cluster formations in the system. As we will see soon below, such phase transition from the gaseous state 
to the cluster state actually appears and characteristic structures are realized in this intermediate energy 
range. 



B. Three phases in SGR model 



As is seen in Fig.H, there are apparently three phases according to the energy per particle U; a) low 
energy collapsed phase (C -phase) U < C^c_iow(^)' ^) intermediate energy phase (X -phase) t^c-low(^) "^ 
U < t/(,_wj^(e), and c) high energy gaseous phase {G -phase) t/p_high(^) "^ ^- "^^^ ^ -phase (c) is 
stabilized by the infrared cutoff {9 < 2tt, or a largest physical scale ^ R), without which the particles 
would escape into spatial infinity. The C -phase (a) is stabilized by the ultraviolet cutoff e, without which 
the particles would fall into a singularity. The specific heat for these particles (a and c) is positive, in 
accordance with the stability of these phases. On the other hand, the I -phase (b) (the most specific to 
gravity) has negative specific heat and therefore is unstable. The nature of this phase is independent of 
any artificial cutoffs and therefore is thought to represent intrinsic properties of gravity. 

The existence of three different phases, as is explained above, distinguishes the SGR model from other 
models with long-range force. For example, HMF model has only high and low energy phases, and shows 



a second-order phase transition between these phases. On the other hand in SGR model, the intermediate 
phase with negative specific heat exists and is strongly unstable when the system is in contact with a 
heat bath. As for IDS model, only a single phase exists because the system has no characteristic energy 
scale; there is no phase transition in this model. 

The limit e — > in SGR model represents genuine gravity without collision and the range of the 
intermediate phase increases without bound. 
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FIG. 3. The time evolution of 7(r) defined by the left side of (A5) in the case with U = -100, -300, -500 



and e = 2.5 x 10 .In each case, the system achieves virial equilibrium within a few dynamical time r. 



C. Three species: gas , halo , and core particles 

In each phase, the particles of the system prevail in various energy ranges. For example, in the 
intermediate phase, some particles evaporate from a cluster and move along the ring almost freely with 
the time scale ta, while some of the others fall into the center of a cluster and oscillate with the time 
scale ij. Thus, the overall phase information does not precisely specify the nature of individual particles. 
In order to obtain much fine information, we define three species in each phase by using the energy of 
the particles: The energy of the I'th particle is given by 
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^ V2iV^l-cos% +e' 



(3.4) 



The classification is based on particle energy, that is, a) core particles for Ei < C/„_1q (e), b) halo particles 
< Ei < C^c_liigh(^)' ^^'^ *") 9'^^ particles for t/c-high(^) < ^i- U 



for t^c-low(«= 



** Note that even in the C -phase, there exist a few gas particles. We term the core particles and halo particles 
as a cluster in this paper; they form an apparent single bound state. 
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FIG. 4. The time evolution of the percentage of each species; Xcore, ^fialo' ^^"^ ^gas in high-energy {Q -) 
phase (top), in intermediate {X -) phase (middle), and in low-energy (C -) phase (bottom) for e — 2.5 x 10~^. Gas 
particles dominate in the high-energy phase, and core particles dominate in the low-energy phase. On the other 
hand, all three species coexist in the intermediate phase. 



In the low-temperature phase of SGR, the quasi-equilibrium state at very low temperature is highly 
inhomogeneous and a single cluster is formed. Most particles are condensed in this cluster and the total 
potential is deep. In the intermediate phase, many particles spread around a cluster. As the temperature 
increases within this phase, halo particles gradually dominate core particles and eventually there appear 
gas particles, which evaporate from a cluster and go along the ring. In the high-temperature phase, all 
particles move almost freely without forming any cluster. 

and Xgas, with A core + 
^. From this, we observe 

0' ^halo ~ 0-2> and 
0.36 for U = -100 {I 



A 



We describe the ratio of particle number in each states Fj by A 



core: 



A 



halo' 



halo 



+ A 



gas 



1. Their evolution is shown in Fig.^ for the case e = 2.5 x 10 



those ratios seem to have approached asymptotic values beyond r w 1: A core 
A^as « 0.8 for U = 100 {G -phase), Acore « 0.41, Xfi^lo ~ 0.23, and Xgas 



" We term a state to represent the situation that a particle belongs to one of the three species. 



-phase), Acore ~ 0.96, Xf^Q^ig ~ 0.02, and Xgas ~ 0.02 for U = —500 (C -phase). Gas particles dominate 
in Q -phase and core particles dominate in C -phase in number. While in I -phase, all of the three species 
of particles coexist almost equally. According to our various calculations changing the particle number 
and initial conditions, this coexistence seems to represent the prominent property of SGS and not the 
finiteness of particles. 
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FIG. 5. 



The middle diagram shows the energy dcpciiucni^c 

-7 



ui uiic iduiu ui each states at time t — 5 for the 
cutoff e = 2.5 X 10~'. For three typical cases, U = 100, —100, and —500, the ratio of time duration intervals 
of the three states in which each particle stays until r = 5 is shown in the upper and lower diagrams. In these 
diagrams, each vertical bar corresponds to a particle, and each ratio of gas , halo , and core state is represented 
by white, gray, and black area, respectively. In both Q -phase {U = 100) and C -phase ((/ = —500), most particles 
stay within each state. On the other hand, in the I -phase {U — —100), many particles experience at least two 
states. Some particles wander in all three states. 



In Fig.^, the energy dependence of the ratio of each species in t = 5 for e = 2.5 x 10^^ is shown at the 
center. The energy scale where halo particles exit corresponds to the one where negative specific heat 
appears (X -phase). Together with it, the relative ratio of time intervals of three states in which each 
particle stays until r = 5 is shown. As we expect, in both Q - and C -phase, most particles stay just in 
one state for quite a long time. For the X -phase, however, many particles experience at least two states. 



Some particles wander from one state to another in three states. 

IV. PARTICLE MOTIONS 

In this section, we examine individual particle motion and velocity distribution function in each phase. 

A. Recurrent motion of halo particles 

In the HMF model at the state near the critical energy, Levy-type flight and anomalous diffusion of 
particles have been reported |13]. These peculiar behaviors of particles are apparently caused by the 
transitions from core particles to gas particles and vice versa. We observe that such peculiar behaviors 
in HMF model are caused by the periodicity of the configuration space and not by the long range force 
itself. In fact, for all Levy- type flights numerically shown in HMF, the flight distances turn out to be 
longer than the period 27r. Therefore the artificial periodicity in the potential, and not the long-rage 
nature of the force itself, is thought to have caused the Levy-type flight in HMF model. In our SGR 
model, we do have the possibility to observe the same Levy-type flight motions since particles in SGR 
also move along a closed ring. Since we would like to extract the intrinsic property of the long-range 
nature of the gravity itself, we pay attention to the recurrent motion of halo particles and disregard the 
round-trip motion along the entire ring. 

First, we choose the parameters as ^ ~ 7.1 x 10~^i? and U — —100 for which the system is in the 
intermediate energy phase. The particle motion is shown in Fig.O(a), in which core particles form a firm 
cluster and they oscillate around the center of the cluster with the time scale t^ and halo particles go in 
and out of the core region without any typical time scale and amplitude. Zoomed in ten-times (Fig.O(b)), 
and even in hundred-times (Fig|2(c)), the particle motion is always similar recurrent movements. These 
repetition of the similar recurrent pattern suggests a self-similar structure of the system. 

Note that the recurrent motion of halo particles is confined within the range [0, 27r] and particles never 
experience a round-trip along the ring, quite contrary to the HMF model. Moreover, this recurrent 
motion of particles is quite robust and is observed in any region of the intermediate energy phase. This 
robustness is a remarkable contrast to the HMF model in which such motion is observed only at the 
critical point in the phase diagram. F^ 

In order to analyze this behavior more quantitatively, we examine the frequency distribution of the 
recurrent time Tj-ec = Tin — ''out, that is, the time period from the moment Tout when a particle leaves 
the barycenter of the core to the moment Tjn when it first returns to the barycenter again. Here we have 
defined the location of the barycenter of the core as 

We 

where 9i is the location of i'th particle in the core and Nc is the total number of core particles at each 
moment. 



■'■'■ The smaller the cutoff e, the larger the range of the total energy per particle U where the recurrent motion 
appears. 
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FIG. 6. Trajectories of halo particles for the time interval (a) r = 5 ~ 8, (b) r = 5 ~ 5.3 and (c) t = 5 ~ 5.03. 
Ten-times zooming up the marked square region in (a) yields (b). Further ten-times zooming up the marked 
square region in (b) yields (c). The recursion profiles are similar despite the scale difference. 



In Fig.^ we depicted the frequency distribution of the recurrent time Tree for core particles and halo 
particles separately. For core particles, as we expected, we find almost Gaussian distribution around the 
center r^ (Fig.0(a)). On the other hand for halo particles, although a peak is found around te{= tE/tji), 
the distribution shows that a long tail spreads widely characterized by the power law /(rrec) ~ '^rec^ with 
p « 2.0 (Fig|j(b)). We have checked that this power law range —3 < logrrec < —1.5 is consistent with 
the eye-fitted region of self-similar motion in Fig.ra. We will see that the above difference of frequency 
distribution for core and halo particles leads to the difference in relaxation time of them in Secjv|. 

We have also examined the energy dependence of the above power p. The result is given in FigJj(b) 
showing that the value of p (« 2.0) is almost independent of the choice of the energy U throughout the 
intermediate energy phase. 
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FIG. 7. A linear-log plot of frequency distribution /(rrec) for core particles and a logarithmic plot of frequency 
distribution /(rrec) for /la/o particles in several values of U. Each plot in (b) shows the same slope within the 
range lO""^" < rrec < 10^^'^. The broken line represents a line with slope —2.0. 



B. Velocity distribution of particles 

As is shown in Appendix |A|, the time averaged kinetic energy {K) is, in general, expressed as the sum 
of iV^ independent stochastic variables. However, in the case of core particles, which are located within 
a few cutoff scale, the quantity (K) turns out to be expressed as the sum of N independent stochastic 
variables (see Appendix H). In this section, we show that the A^-dependence of (K) characterizes the 
velocity distribution of particles. 

In Fig.@, we show the velocity distributions of core and halo particles in three cases; (a) ^ ss 5 x lO^^R, 
U = -100 (b) ^ « 5 X lO-^i?, U = -500 and (c) C « 5 x lO'^i?, U = -0.65. We have superposed 
the velocity data at r = 1, 2, 3, 4, and 5, with ten different random initial conditions fixing the total 
energy. Thus the size of the whole data we used is 5 x 10^ for each velocity distribution function. The 
distribution of core particles is well fitted by the Gaussian distribution 



P{v) 



1 



2'Ka 



(4.2) 



with the dispersion (v^l^^^"^ = ct = 7.29(a), 10.9(6), and 0.51(c). 

On the other hand, the distribution of halo particles in (a) is manifestly non-Gaussian, although that 
of the same particles in the case (c) is Gaussian. This distribution is at least unchanged until t = 200 . 
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FIG. 8. The linear-log plot of velocity distributions of particles for the three cases; (a)e = 2.5 x 10 
{I -phase), (b)e = 2.5 x lO"'', U = -500 (C -phase), and (c)e = 2.5 x 10-^ U = -0.65 (I -phase), 
velocities of core particles are well fitted by the Gaussian distribution with dispersion a — 0.51(a), 7.29(b), and 
0.51(c) respectively. Also in the large cutoff case (c), velocities of halo particles are well fitted by the Gaussian 
distribution a = 0.93. On the other hand, in the small cutoff case (a), velocities of halo particles are well fitted 
by the Lorentzian distribution (q = 1) with s — 2.2. 



What is the origin of the above non-Gaussian distribution for halo particles in the intermediate energy 
phase in case (a)? Here we analyze this issue from the viewpoint of the generalized central-limit theorem. 
It is well known that the limit distribution for sums of independent random variables is the Gaussian 
distribution provided that the dispersion is finite (the central- limit theorem). However, it is less well- 
known that the limit distribution for sums of independent random variables is the stable distribution in 
general cases, including those where the dispersion is divergent. This stable distribution is defined to 
satisfy the relation 



N 



TVV" 



X, 



(4.3) 



where — means that the distributions of both sides are equal to each other. The parameter a classifies 
the stable distributions and must satisfy < a < 2 for the normalizability and the positivity of the 
probability distribution function. In the definition above, Xi and x are the probabilistic variables obeying 
the same distribution. Let us apply this to the velocity distribution of our model. We square both sides 
of the above equation, and obtain the non-extensive property. 



\ /A'' — particle 



\ / one — 



one — particle 



(4.4) 



and the method of characteristic function yields the explicit form of the distribution function ||20|| 

I f°° 
Piv) = 7r dy eM-ivy ~ s\y\'^]. (4.5) 



Note that this stable distribution includes the Gaussian distribution (4.2) as a special case a — 2, where 
the kinetic energy becomes extensive, and the dispersion a = \/2s is finite. The parameter s is thought 
to be a generalized temperature. 

Let us first consider the velocity distributions of core particles. We observe, from our numerical calcu- 
lations, that all core particles oscillate within the narrow region of the cutoff size ^, and the gravitational 
two-body interaction is dominated by the artificial potential force. Therefore the averaged kinetic energy 
is described as the sum of N independent statistical elements (Appendix pi) . This leads to the normal ex- 
tensivity for the velocity distributions of core particles and therefore we expect the Gaussian distribution. 
This is consistent with the results in Fig.R 

On the other hand, for halo particles which interact through genui ne g ravity, the kinetic energy behaves 
as (Kp) ^ N^ for fixed R as is shown in Appendix ^. Thus from (44) we fin d th e index a = 1 for the 



physical velocity distribution of SGR. In this case with a — 1, the distribution ( [4.5D becomes the Lorentz 
form 

Pi^)--^^- (4-6) 

This is also consistent with the results in Fig.^. The essence of the appearance of this non-Gaussian 
distribution is the non-extensivity of the energy for SGS; this is the intrinsic property of gravity. 

Then what is the origin of the Gaussian distribution of halo particles in case (c)? As we show in the 
next section, this is mainly because the relaxation times of core and halo particles are close with each 
other. In the next section, we will discuss the relation between the profile of velocity distribution and 
the relaxation process. 



V. RELAXATION 

One of the most important issues in statistical physics of the A^-particle system is the relaxation process. 
In order to study the relaxation process of our model, we choose several different initial conditions for the 
same values of e and U. Here we have examined two parameter sets, i.e. (i) e — 2.5 x 10""^, U — —0.65 
and (ii) e — 2.5 x 10^^, U ~ —100, for both of which negative specific heat appears. 

As we discussed in SecO, there are three dynamical time scales: tn, t^, and tE- In the units of our 
normalization, the ratios of these time scales are 

tc/t„ = {2tf'\ (5.1) 

For example, if we choose U — —100 and e — 2.5 x 10^^, then tE/ta ~ 1.8 x 10^^ and therefore the three 
time scales are separated as 

h<tE<: tR. (5.2) 

Then the contribution of gas particles to the dynamics of halo or core particles is negligible during a few 

tR. 



The parameter set (i) does not meet the condition ( |3.3|) and IeI^r ~ 0.34; Ie is the same order as Ir 



On the other hand, the parameter set (ii) meets (3.3) and tEJtR ~ 1-8 x 10 , Ie is much shorter than 



Ir. This means that gas species play an important role in relaxation in case (i), while they do not in case 
(ii). In order to examine the relaxation process for these two cases, it is convenient to normalize the time 
with the unit t^, since the degree of relaxation of 3-D gravitational system is conventionally measured 
with the dynamical time. So here we introduce the dimensionless time te = t/tE- 

As for initial conditions, we locate several clusters of the same size with the same interval. The number 
of clusters is chosen as 5, 10, 20, 25, and 50 for the model (i) and 5, 10, 20, and 25 for the model (ii). 
The initial velocities of particles are set randomly but the total energy is fixed: U = —0.65 for (i) and 
[/ = -100 for (ii). 

In order to study the Ergode property of the system, we examine the mixing property of the particles 
in the three states( core , halo , and gas ). In other words, we examine the degree of isolation of each 
states. We directly measure how extent each particle experiences these three states. First, we pick up 



discrete times {tei, te2, • • • , tei, ■ ■ •) with the equal interval Ate', at each time, each particle stays at one 
of the three states. Then we count how many times the /c-'th particle stays at each states before the time 
TEi- Finally, after normalization, we obtain the relative frequencies of three states AgQy.g(rEi), A^ , {'TEi) 

and Xgasi'^Ei) (with Xcorei'^Ei) + -^/ia/o('^£;j)+A|as('rsj) = 1 ) for the fc'th particle. 

At the first time tei, each particle definitely stays at one of the three states; relative frequencies of 
stay Ago^g(T£;i), A^ , (tei) and X'gQ^glTEi) are either 1 or 0. In the next time te2, some particles may 

change the state; then X'^orei'''E2) , X'f , {te2) and XgasiTE2) for these particles are 1/2. In this way, 
^corei'^Et), X'l^lgiTEi) and X^gasi^Ei) wiU evolve in time t^j. Various distribution of {X'l.oreiTEi)}k=i 
for all the particles k = 1,2, ■ ■ ■ , N will define the distribution function N{Xgas, te) for values Xgas 
at time tei by simply counting the number of the particles which take the value Xcore at time te- 
Similarly distribution functions for other states can also be defined: N{Xi,„i„, te) and N(Xcore, te)- 
As the system is thermally relaxed and the equipartition of energy is attained, AgQj~g(r^), A^ , {te) and 
Xqasi'^E) will converge to the same value Ag^j-g, At , and A^^^ independent of the particle label k. Then 
N{Xgas, Te), N{Xf^^ig, te) and N{Xcore, te) would have sharp peaks at X*coreiX\^i^ and A^^^ in the 
thermal equilibrium state. Thus the development of the peaks in the distributions functions N[Xgas, te), 
N{Xf^^lg, Te) and N{Xcore, te) can be a good measure of the degree of thermal equilibrium. 

In the large UV-cutoff case (i), each distribution N{Xgas, te), -^(A/j^/o' '''^) ^^'^ N{Xcore, te) starting 
from Tei — 6.5 x 10* seems to develop a single peak in-between and 1 at least by the time te ~ 10^ 
(Fig.y (b)), and the variance around the peak seems to reduce in time (Fig.|| (c)). This means that 
most particles experience all of the three states at least by te ~ 2.0 x 10^; thermal relaxation proceeds 
in the time scale te ^ 10^. Contrary in the small UV-cutoff case (ii), each distribution N{Xgas, te), 
N{Xf^^lg, Te) and N{Xcore, te) starting from tei = 1.1 x 10^ seems to develop peaks at the edge of 
the domain 1 and/or 0, even after the time duration, te ~ 3.8 x 10^(Fig|9| (e)) . This means that most 
particles stay in each single state even in the time scale te ~ 10^; the system does not approach thermal 
relaxation at all in this time scale. 
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FIG. 9. The time variation of A'^(A, te), where the bar in white, gray and black colors show N{\gas, te), 
^(■^halo' '''^) ^^'^ N{Xcore, te), respectively. The left and right columns respectively represent the case of large 
UV-cutoff (a) e = 2.5 X lO-^* and U = -0.65 and the case of small UV-cutoff (b) e = 2.5 x 10"^ and U 



where tei = 6.5 x 10'', Ate = 2.0 x 10^ for (a), and tei = 1.1 x 10^, 
time flows from top to bottom: (al, fol) tei , (a2, 62) te^ , (a3, 63) tei2 i 



-100, 
Ate = 3.8 x 10^ for (b) In both columns, 
(a4, 64) TE34 , respectively. 



By combining these results, the Gaussian velocity distribution reahzed in the case (i) in Fig.H seems 
to reflect the achieved thermal relaxation of the system. On the other hand, the Lorentzian velocity 
distribution realized in the case (ii) in Fig.g seems to reflect some quasi-equilibrium intermediate stage 
before the thermal relaxation of the system. Latter peculiar properties in velocity distributions may make 
3-D gravitational system quite different from those in the system with short-range interactions or the 
system with positive specific heat. If we set e larger, then the encounter of halo particles with normal 
particles {gas and core particles) becomes more frequent, and the halo particles cannot be isolated from 
the normal particles, resulting trivial thermodynamical equilibrium of the whole system as in the case 

(i)- 



VI. SCALING PROPERTIES OF THE HALO 



In Sec.|rv|, we have found that the velocity distribution of halo particles is non-Gaussian with anomalous 
w~^-tail and each halo particle shows intermittent recursive motion around the core without a definite 
time scale. These scale-free properties of the halo particles turn out to appear in the spatial distribution 
profile of them itself. In this section, we study the scaling property of the spatial distribution with the 
box-counting method pl| . 

For this purpose, we use the box counting method: We divide the entire configuration space 2t:R into 
segments with equal size i and count the number N{£) of the segment which contains at least one halo 



particle. Then we define the quantity — 9(log A^(^))/9(log^) which turns out to be the scahng exponent, 
provided that the quantity is almost independent of the scale L In order to increase the statistical 
significance, we used multiple particle-distribution data at difi:erent times for each run. Moreover, we 
made many runs of calculations with different initial conditions for the fixed total energy U . We have 
applied the box-counting method to each data and then all the individual results are superposed. When 
we extract the distribution data, we have chosen the time interval [2000t^, 5000t^], in which the particle 



ratios Xt 



A 



halo^ 



and A gas almost reach the relaxed constant values. Thus we calculate the following 



averaged quantity 



D{i) 



d{\ogN{e)) \ 

d{\ogl) /^ 



(6.1) 



over all the data we thus prepared. If any scaling property exists in the particle distribution, D{£) would 
become constant for a finite range oi £i < £ < £2- As for the bounds £1 and £2, we have technical 
restrictions originating from our numerical calculation method. The averaged distance between gas 
particles is estimated as 



%as(e) 



2TrR 



NX 



gas 



(6.2) 



If the box size £ is larger than ^gas(e), we expect D ^ 1.0, because almost all the boxes contain at least 
one particle. On the other hand, if £ is smaller than the cutoff distance ^, the genuine property of gravity 
is lost. Therefore scaling property in particle distribution is relevant only within the rage ^ < £ < Igas- 
Actually in our calculation, for e = 2.5 x 10~^ and U — —100, we can see the typical behavior of N{£) as a 
function of £ in Fig.|l^. In this plot, we can see scaling behavior N{£) ~ £~'^ with the small exponent in the 
range lO^'^i? < £ < 10~^R, which is well inside the relevant region. This scaling seems to originate from 
halo particles, because they not only dominate in the above scaling range but also show non-Gaussian 
velocity distributions and self-similar recursion jumps. 

However, the scaling property of the box counting method does not always conclude the existence of 
the fractal structure in the particle distributions. This is because the box counting method itself cannot 
distinguish the two possibilities; a) genuine fractal structure and b) the power-law tail of the particle 
distributions around the core center. 

In order to distinguish the above possibilities, we compare two different superposition methods as 
follows: a) We simply superpose all the data with the bare coordinate 9 and b) we superpose all the data 
with the coordinate adjusting so that the mass-center of each data comes to the same position; i.e. we 
introduce a new coordinate 9 as 9 = 9 — 9bc for each data with the mass center 9bc- 

As a result, the scaling property of the particle distributions, which was observed in the original 
method, disappears completely in the case a). On the other hand in the case b) we observe the same 
scaling property as was observed in the original method. This result suggests that the observed scaling 
property in particle distribution is not the genuine fractal property of the system, but simply due to 
the power-law distribution of particles around the mass-center. We have confirmed this with our further 
analysis on the direct observation of the power-law particle distributions around the mass-center. 
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FIG. 10. Log-log plot for the box scale L vs the number of the box A'' occupied by at least one halo particle. 
We superposed the snapshot data of n different times. The cases with n = 1, 10 (to = 2.0, Ar — 0.3), and 



100 (to = 2.0, Ar = 0.03) with the cutoff e = 2.5 x 10"' 
cases coincide with each other in the range ^ < L < 10 
in this range is 0.1. 






7.1 X 10 R) are plotted. The plots of these three 
The scaling exponent D{L) derived from the slope 



The energy and cutoff dependence on the scaling exponent D is shown in Fig.^. Here we show the 
scaling exponent of halo particles in the region of negative specific heat. We find that the scaling exponent 
D is almost constant in this region. Moreover, smaller cutoff e corresponds to the smaller dimension D. 
On the other hand in the limit of e ^ 1, the scaling property cannot be observed. 

In the limit of real gravitational interaction (e —^ 0), the region of negative specific heat extends in 
energy range and so does the scaling region. This constant property of the exponent D is likely to be 
independent of the temperature or the energy of the system. 
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FIG. 11. The scaling exponent D vs the energy per particle U. Most of the energy region shown in the figure 
corresponds to T -phase. For convenience, we use y/eU as a horizontal axis. The case of e = 2.5 x 10~^ (squares 
with error bars ) and the case of e = 2.5 x 10~^ (crosses with error bars ) are plotted. It seems that the dimension 
D does not depend on U. Since the number of halo particles is small in low energy, the error becomes large there. 
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VII. CONCLUSION 



In this paper we have studied the self-gravitating ring (SGR) model which is a one-dimensional (1-D) 
model of a self-gravitating system. Although the system is 1-D, the particles constrained on a circular 
ring mutually interact with the 3-D full gravitational force. We found that an interesting phase with 
the negative specific heat appears at the intermediate energy scale, which reflects the virial condition of 
the 3-D gravitating system. Classifying the particles in each phase into three species {core , halo , and 
gas specie), we examined dynamical properties of this phase. Using SGR model, we further examined 
characteristic properties of relaxation, velocity distributions, and density profiles in each phase. 

The cutoff parameter ^ introduced in SGR model controls the singular properties of the system at 
short-distances. By changing this parameter, we could systematically study the effect of the short- 
distance singularity of the force upon the statistical properties of the system. 

When the cutoff parameter is of order one (^ ~ i? ), the system resembles the Hamilton Mean Field 
model (HMF). In fact, the velocity distribution becomes Gaussian at i ~ 0(10'^) t/j. This is because the 
two time scales Ie and t^ become the same order and therefore halo particles frequently encounter and 
interact with gas particles. This strong interaction induces frequent exchange of particles between these 
two species. Thus the energy diffuses effectively and complete thermodynamical equilibrium of the entire 
system is established. Especially the temperature is the same for all the particle species. In other words, 
the negative specific heat region {halo particles) does exist, but the strong interaction with gas and core 
particles (positive specific heat) quickly dissolves the negative specific heat region. 

In much more smaller cutoff case with ^ w 7.1xl0"^i?(e = 2.5x10"'^), the velocity distribution becomes 
almost Gaussian after the virialization, although the temperature of each species is different from each 



other. This temperature difference, as well as the required time for establishing the equilibrium, is 
increased for smaller value of the cutoff £,. 

On the other hand, as ^ ^ 0, the difference of the two time scales increases Ir ^ t^, and the energy 
region of the negative specific heat extends wider. Therefore we can expect that halo particles become 
almost isolated from the gas particles. This isolation of unstable halo particles from other normal species 
enables the halo species to last sufficiently long. In the case of ^ w 7.1 x 10~^i?(e — 2.5 x 10~^), the initial 
condition dependence of halo particles survives even after 5 x IOHr in our simulation, which prevents 
the full system from reaching thermal equilibrium. These results resemble the case of core particles 
with <^ « 7.1 x 10^^ R. However, the remarkable difference between the two cases is that the velocity 
distribution of halo particles in ^ « 7.1 x 10~^i? is not Gaussian but Lorentzian distribution, while that 
of core particles in ^ « 7.1 x 10~*i? is Gaussian. This property of halo particles qualitatively reflects 
genuine 3-D gravitational interaction and is the main difference from the HMF models or IDS models We 
have also found the scaling structure of halo particles in the intermediate energy phase. When we reduce 
the cutoff ^ ^- 0, the scaling region as well as the energy range of negative specific heat is extended wider, 
although the value of exponent depends on the cutoff and reduces with ^. 

In IDS model, the system reaches thermodynamical equilibrium at least at t ~ 0(10^)^^;, where i^ = 
(l/47rGmiV) ^jAE/mN is the crossing time of IDS model [Q. It is also shown that the system reaches 
the collisionless mixing phase during much shorter time interval. Our simulation with ^ « 7.1 x IQ^^R 
shows non-Gaussian velocity and scaling property in such a short time interval 4 x 10^ t^; ~ 2 x 10^ Ie for 
U = —100. So the non-Gaussian velocity distribution might reflect the character in collisionless mixing 
phase which appears in IDS simulation |11[| . However, it might also be true that the relaxation process in 
SGR is quite different from that in IDS, where no time scale separation between core and halo particles 
exists. So the exotic character of halo particles might be intrinsic in SGR or the system with real 3-D 
gravitational interactions because of the time scale separation. 

In summary, SGR model is characterized (i) by the phase separation and the particle state separation 
and (ii) by similar recurrent motion and non-Gaussian velocity distribution of halo particles . Both of 
these characters may play a key role in explaining the observational properties in real 3-D self-gravitating 
objects. For example, the recent observations with Hubble Space Telescope support the existence of 
supermassive black holes or cusps at the center of triaxial elliptical galaxies [p3[-p5[. This observation 
suggests that the above property (i) may appear in these galaxies, since they have high density region 
corresponding to core particles at their centers and low density region corresponding to halo particles 
surrounding them. Hence, the property (i) may give a hint for explaining the relaxation process and 
the origin of stationally conflguration shape of these galaxies, since the stochastic halo orbits which 
interact with core particles at the center of a triaxial galaxy may affect the equilibrium shape through 
the continued mixing near the center |2q] . In addition, the propety (ii) might help to explain the observed 
fractal structures and non-Gaussian velocity distributions in the inter stellar medium p7| , which are likely 
to be gravitationally virialized. 
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APPENDIX A: VIRIAL CONDITION FOR SGR MODEL 

We start with the non-dimensional Hamiltonian of the SGR model, 

N N 



2tt ^%/2A^v/l-cos%+6 



Canonical equations of motion for this system become 



dp. 






dT j^^2^/2N{l-cos^^J+e)^/^' 



Then the second time derivative of the inertia moment of this system becomes 

i—l \^— 1 2—1 / 



^/' f^ V2A^(1- COS %+ 6)3/2- 
Taking long-time average for this equation, we obtain 

\dr'tl / \^ ^/27V(l - cos % +6)3/2/ 

Therefore if the condition 

/ r/2 ^ \ 



dr2 

I 4=1 



is reahzed, then the virial relation 



2iV \^ V2N{1 - cos e,, + e)3/2 / ' 

holds. 

By using the total potential 

^ 1 

y^-Y^ , (A7) 

f^ V2N^l-cos9,,+e 



the virial relation (A6) can be expressed as 



2N \'^ de, 



Note that even in the limit of e ^ 0, the virial relation {K) = 2NT = — {V) /2 is not justified, because 
of the potential form (A7). However, as is shown in SecO, we can get the above ideal virial condition 



for 3-D gravity within the appropriate energy range. 

APPENDIX B: VIRIAL CONDITION IN A C -PHASE FOR SGR MODEL 

In C -phase, almost all of the particles are trapped inside the core , whose size is almost ^. Therefore 



in this state canonical equations (A2) are approximated as 

de, 
17 ^P" 



Jtj 



dpi _ Y^ 



1 \ ^ f) 

1-^ +E:r7&7^> (Bi) 



2V2e3/2 V N ^. 2V27Ve3/2 



since |%| ^ 0{e) < 1. 

In this approximation, the second time derivative of the inertia moment of this system becomes 

i i ^ ^ '^ i—1 ^ i<Cj 

Taking long-time average for this relation, we obtain 



Since 6i and 9j{i ^ j) are independent stochastic variables, the last term of (P33| ) vanishes. Hence from 
the virial relation, 



T 



and 



^-^)tm^ m 



2V2e3/2i? V N 



i=l 



N 



('«-:;:Sk; '4 EC?). C5) 



i=l 



hold. This shows that (Kp) in C -phase is given by the sum of N statistical elements. 
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